Non mean-field behavior of the contact process on scale-free networks 
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We present an analysis of the classical contact process on scale-free networks. A mean-field 
study, both for finite and infinite network sizes, yields an absorbing-state phase transition at a 
finite critical value of the control parameter, characterized by a set of exponents depending on the 
network structure. Since finite size effects are large and the infinite network limit cannot be reached 
in practice, a numerical study of the transition requires the application of finite size scaling theory. 
Contrary to other critical phenomena studied previously, the contact process in scale-free networks 
exhibits a non-trivial critical behavior that cannot be quantitatively accounted for by mean-field 
theory. 
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Recent years have witnessed an upsurge in the study 
of the structure and function of complex networks 
The interest in these systems has been driven by two 
motivations. First, a large amount of work has been de- 
voted to the structural characterization of real networks, 
the main contribution being the discovery of their seem- 
ingly ubiquitous scale- free (SF) nature Q, characterized 
by a probability P(k) that an element in the network 
(a vertex) is connected to other k elements (has degree 
k), scaling as a power-law, P(k) ~ fc~ 7 , with a certain 
degree exponent 7 usually in the range 2 < 7 < 3 Q. 
On the other hand, a great deal of attention has been 
drawn by the profound effects that a heterogeneous SF 
connectivity pattern can have on dynamical processes on 
complex networks. In particular, models exhibiting criti- 
cal phase transitions such as ferromagnets 0] , percolation 
or epidemic models p| have been studied on this class 
of networks. Present understanding of critical behavior 
in complex networks is essentially based on the applica- 
tion of mean-field (MF) techniques [||. In the study of 
critical phenomena on regular substrates (lattices), MF 
theory plays a prominent role. Despite its simplicity, 
it describes qualitatively most phase transitions. More- 
over, while in low spatial dimensions fluctuations cannot 
be neglected, above the upper critical dimension MF is 
even quantitatively correct, providing exact expressions 
for critical exponents. In the case of complex networks, 
since they usually fulfill the so-called small-world (SW) 
property 0, such systems can be viewed as infinite di- 
mensional, and therefore one can expect MF to provide 
an exact description of dynamical processes taking place 
on networks, once degree heterogeneity has been taken 
into account. Results obtained so far seem to validate 
this idea, particularly for SF networks [8| . 

In this paper we show the inadequacy of MF theory to 
describe some dynamical processes on SF networks, by 
considering the absorbing-state phase transition of the 
contact process (CP) 0. We develop a MF theory for 



the CP on generic networks, and solve it for SF net- 
works, obtaining expressions for both infinite and finite 
networks. These expressions are compared with numer- 
ical simulations of the CP on SF networks, which are 
analyzed applying a finite size scaling (FSS) theory [l^. 
Our results show that the the critical exponents of CP 
on SF networks cannot be accounted for by MF theory, 
and signal the presence of nontrivial fluctuation effects. 

We consider the classical CP on generic networks, de- 
fined as follows : An initial fraction po of vertices is 
randomly chosen and occupied by a particle. The time 
evolution of the process runs as follows: At each time 
step t, a particle in a vertex Vi is chosen at random. 
With probability p, this particle disappears. With prob- 
ability I— p, on the other hand, the particle may generate 
an offspring. To do so, a vertex Vj, nearest neighbor of 
the vertex Vi, is randomly chosen. If vertex Vj is oc- 
cupied by a particle, nothing happens; if it is empty, a 
new particle is created on Vj. In any case, time is up- 
dated as t — » t + l/n(t), where n(t) is the number of 
particles present at the beginning of the time step [Io| . 
In Euclidean d-dimensional lattices, the CP undergoes a 
continuous transition at a critical point p c , separating an 
absorbing phase from an active one 0] . This phase tran- 
sition is characterized in terms of the order parameter p, 
defined as the average density of particles in the steady 
state. For p > p c , an absorbing phase with p = is ob- 
served. For p < p c the system sets in an active phase with 
a nonzero order parameter, obeying p ~ A' 3 for p — > p~ , 
where A — p c ~p. Close to p c , the system is further char- 
acterized by diverging correlation length and time scales, 
namely £ ~ lAj^- 1 - and r ~ |A| _I/ ii. The critical expo- 
nents (3, v± and v\\ provide a full characterization of the 
phase transition of CP in Euclidean lattices. Below the 
critical dimension d c = 4, the exponents are nontrivial, 
and depend explicitly on d. For d > d Cl the exponents 
take the classical MF values (3 = v\\ = 1, v±_ = 1/2. 

In order to gain analytical insight on the behavior of 
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the CP on complex networks, we consider as usual a MF 
rate equation for thepartial densities Pk(t) of occupied 
vertices of degree k (5|, from which the total density is 
obtained as p(t) = ^2 k p k (t)P(k). The rate equation 
for the partial densities in a network characterized by 
a degree distribution P(k) and degree correlations given 
by the conditional probability P(k'\k) that a vertex of 
degree k is connected to a vertex of degree k' 0] can be 
thus written as 0] 



d tP k(t) = -p k {t) + Afc[l - p k {t)] 



P(k'\k) Pk ,(t) 



(1) 



The destruction term, due to the vanishing of particles — 
with probability p — at vertices of degree fc, is propor- 
tional to the density of occupied vertices of degree k, 
p k (t). The creation term, on the other hand, is due to the 
generation of offsprings from particles located at nearest 
neighbor vertices of empty vertices of degree k. Thus, 
it is proportional to the probability [1 — p k (t)] that a 
vertex of degree k is empty, that it is connected to a 
vertex of degree k' , P(k'\k), that this one is occupied, 
pk'{t), that it generates an offspring, 1 — p, and that this 
neighbor chooses to create the offspring on the vertex of 
degree k under consideration, (1/fc'). In Eq. Q we have 
performed a rescaling of time, introducing the new pa- 
rameter A = (1 —p)/p. For uncorrelated networks, where 
P(k'\k) = k' P(k') I ' (k) [llj, the rate equations take the 
form 



d tPk {t) = -p k (t) + X—[l 



p k (t)]p(t). (2) 



The solution of Eq. (J2J) depends on the nature of 
the degree distribution. For homogeneous networks, in 
which P(k) decays exponentially for large degrees, all ver- 
tices are approximately equivalent. Therefore, we have 
k ~ (fc) and p k (t) ~ p(t), so that the rate equation reads 
d t p(t) = -Xp 2 (t) + Ap(t), where A = A — 1. This equa- 
tion describes the behavior of the CP in an infinite di- 
mensional system; the critical point A c = 1 (p c = 1/2 ) 
and the exponent (3 = 1 are thus immediately found |10| . 
Moreover, from the solution p(t) ~ A/A — e~ At , we can 
identify the relevant characteristic time r ~ A -1 and the 
corresponding exponent i/u =1. 

For heterogeneous networks, with a degree distribution 
exhibiting large fluctuations, the solution of Eq. Q must 
consider explicitly the k dependence of the partial den- 
sities. Information on the active phase is obtained by 
imposing the steady state condition, dtPk{t) = 0, which 
yields the nonzero solutions 



Pk 



\kp/(k) 
l + Xkp/(k}' 



(3) 



where p k is now independent of time. In order to de- 
termine the critical behavior of the average density p as 



a function of the control parameter A we focus on the 
self-consistent equation for the order parameter p. By 
combining Eq. © with the expression for p, one obtains 



\p X - fcP(fc) 
(fc) V 1 + XkpUk)' 

k 



(4) 



This equation depends on the full degree distribution. In 
the case of SF networks, for which the degree distribu- 
tion in the continuous degree approximation is given by 
P(k) = (7 — l)77i 7_1 /c~' 7 , with m the minimum degree in 
the network, the solution will depend on the degree ex- 
ponent 7. Substituting the summation by an integral in 
Eq. we obtain in the infinite network size limit (i.e. 
when the degree belongs to the range [to, 00]) the expres- 
sion p = F[l, 7—1,7, — (Xpm/ {fc}) -1 ], where F[a, b, c, z] 
is the Gauss hyper-geometric function To evaluate 
the critical behavior for small p, we invert this expression 
using the asymptotic expansion of the hyper-geometric 
function for z — ► —00, obtaining 0] p(A) ~ (A— l) 1 /^- 2 ) 
for 2 < 7 < 3 and p(X) ~ A - l for 7 > 3. Hence the MF 
solution in infinite networks gives the critical exponent 
(3 MF = 1/(7 - 2) for 7 < 3, and /3 MF = 1 for 7 > 3, 
while the critical threshold is A^ F = 1 {p™ F = 1/2), 
independent of the degree exponent. 

The previous results, corresponding to networks of in- 
finite size, are however strongly affected by finite size 
effects. The difference between an infinite SF network 
and one of finite size N is that the latter has a cut-off or 
maximum degree k c , which is a function of the network 
size and that, for uncorrelated networks, scales as 
k c (N) ~ N 1 / 2 0| . The presence of the cut-off restricts 
the possible values of the degree to the range [m, k c ] . 
For finite networks, Eq. Q yields, in the continuous de- 
gree approximation, p = F[l,~f — 1,7, — (Xpm/ (k))^ 1 ] — 
Oc/to) 1 "^^ - l,7,-(Apfc c /(/c))- 1 ]. To obtain the 
MF critical behavior for a finite network, we evaluate this 
expression in the limit p — ^ for generic pk Cl obtaining 



(k) ( (fc)Afer 2 > 

XkJ \ Ato^- 1 



(5) 



where f(x) is the inverse of the function g(x) = (F[l, 7 — 
1,7,-1/x] + r( 7 )r(2 - 7 )xT- 1 )/a;. For x > 1, f (x) 
diverges as x 1 ^ 1 ^, in agreement with the infinite size 
limit p ~ A^ . For small arguments, the function f(x) 
vanishes for x — > x^ with xq > 0. Therefore, for finite N, 
if A is made smaller than a value proportional to k 2 ~~<, 
then p = 0. This indicates that the effective critical point 
is shifted in the active phase in finite networks and that 
the convergence to the infinite network limit is slow, since 
corrections vanish as iV^ 2 " 7 )/ 2 , i.e., with an exponent 
smaller than 1. 

^From Eq. Q we can ascertain when, in a system of 
N vertices, finite size effects start to play a relevant role. 
It is clear that for x < 1, i.e. A < A x ~ (fc c /m) 2 ' such 
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effects start to appear, so that one expects deviations 
from the decay p ~ A p . In terms of the order pa- 
rameter this implies that finite size effects are negligible 
only if p 3> Px = (k)/(Xk c ). This condition introduces 
serious constraints on the precision of the numerical esti- 
mate of the critical point when performed by measuring 
the dependence of the p on p in the active phase. 

This difficulty can be overcome by using the finite size 
scaling (FSS) technique [HE3- The FSS technique is 
based on the observation that, even below the critical 
point, the density of active sites in surviving runs p a 
reaches a quasi-steady state whose average is a decreasing 
function of the network size TV [Tfi| . The FSS ansatz 
assumes the dependence of p a as a function of A = p c ~P 
and N 

p a (A, N) = /(A/V 1 /^), ( 6 ) 

where the scaling function has the asymptotic behavior 
f(x) — ► x@ for x — > oo, and f(x) — ► const, for x — > 
0. Compatibility with the theoretical behavior given in 
Eq. © for x 3> 1 requires that the exponents assume 
the MF values (3 = 1/(7-2) and v± = 2/(7 - 2) in 
uncorrelated networks with cut-off scaling as TV 1 / 2 . Thus, 
at p c the function p a is expected to decay as 

p a (0,N)~N-V Uj : (7) 

This fact can be used to estimate the critical point, as 
the value of p yielding a power-law behavior for p a as 
a function of N. MF theory predicts in this case an 
exponent /3/Vj_ = 1/2. 

In order to check these results, we have simulated nu- 
mericall y th e CP on the uncorrelated configuration model 
(UCM) [lTj , which generates uncorrelated SF networks 
with arbitrary degree exponent and cut-off k c (N) ~ 
TV 1 / 2 . We consider the particular cases 7 = 2.75, 2.5 and 
2.25. Simulations were performed with at least 100 runs 
of the CP for each realization of the network and no less 
than 100 such realizations for each value of 7 and N. 

We focus first on the determination of the critical point 
and of the exponent ratio (3/v± via Eq. J7J. In Fig.^ wc 
display the density of active sites in surviving runs p a as 
a function of the network size, for 7 = 2.75 and several 
values of the control parameter p. The critical point is 
identified as the value of p such that the decay of p a vs 
N is a pure power-law, which for this particular value 
of 7 corresponds to p c — 0.4215. In Fig. we show the 
decay of p a vs N at criticality, for the three values of 7 
considered. In all cases a very clear power-law behavior 
is identified, yielding the exponent ratios (3/v± = 0.63 
for 7 = 2.75, 0.70 for 7 = 2.50, and 0.76 for 7 = 2.25. 
Notice that these values are in strong disagreement with 
the MF prediction 0/u± = 1/2. 

In principle the exponent (3 should be determined by 
the way the steady state density of particles vanishes as a 
function of p c — p in the limit N — > 00. However, it turns 
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FIG. 1: p a (p,N) as a function of the network size N for the 
CP on UCM networks with 7 = 2.75, for values of p close 
to p c . Solid lines are guides to the eye. Inset: p a (p,N) for 
p = 0.36, showing that the infinite TV limit is not reached even 
far from p c . 




FIG. 2: Scaling of p a {0, N) as a function of the network size 
for the CP on UCM networks for three values of 7. The 
critical points and exponents /3/u± obtained are reported in 
Table Solid lines are least squares fits. 



out that, even far from p c , it is numerically impossible to 
fulfill the condition p 3> p x that guarantees the absence 
of finite size effects. This is apparent in the inset of Fig. ^ 
where one can see that even for p = 0.36 the value of p a 
slightly depends on N even for large sizes. This very slow 
convergence hampers the determination of (3 in the usual 
way. 

The numerical value of (3 and v±_ can be however deter- 
mined via the FSS analysis, by noticing that from Eq. © 
we expect data for the surviving density to collapse onto 
a universal function f(x) when plotting N^^ u± p a (A, N) 
as a function of AN 1 / U± . Using the values of /3/^j_ ob- 
tained from Fig. ^ we try the collapse for various values 
of (3 (and, as a consequence, of The best data col- 
lapses, shown in Fig. [31 are obtained for the values of (3 
and u± reported in Table [I] The scaling function f(x) 
is expected to diverge as x@ for asymptotically large x. 
We measure how the numerical scaling functions go for 
large x in Fig. [21 finding the effective behavior f(x) ~ x@ 
with (3' = 1.43, 1.70, and 1.79. The disagreement with 
the values of (3 providing the best data collapse indicates 
that the values of x that can be reached are not large 
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FIG. 3: Data collapse of N 0/ " ± p a (A,N) as a function of 
AN 1/v± for the CP on UCM networks. Data are for system 
size ranging from N = 10 3 to TV = 10 7 . For clarity the data 
for 7 = 2.25 have been shifted upwards. The values of v± are 
those giving the best data collapses. The straight lines are 
least squares fits to the behavior f(x) — > x 13 , for large x. 

enough to see the truly asymptotic x@ behavior. 

We summarize the results obtained in Table [I] The 
conclusion of this analysis is that the CP in SF net- 
works exhibits an absorbing-state phase transition at a 
non-zero critical point, whose exponents are not correctly 
predicted by MF theory [l^ . This result is surprising: it 
is generally believed that dynamical models on SF net- 
works find in this formalism a correct description, once 
degree heterogeneity has been taken into account Q . The 
results presented here open two clear perspectives for fu- 
ture work. The first is the quest for a theory beyond MF 
for the CP on heterogeneous networks. The second is the 
investigation of the limits of MF theory for generic dy- 
namical processes on heterogeneous topologies. A useful 
indication along these lines is provided by Fig. where 
we plot the ratio of the standard deviation Apk over the 
mean density pk restricted to vertices of degree k, at the 
critical point. It turns out that the fluctuations are much 
larger than the average except for the smallest k, all the 
more so for increasing system size. Hence it is incorrect 
to assume (as MF theory does) that the mean pk pro- 
vide all information about the distribution of densities 
for vertices of a certain degree. The inclusion of these 
large fluctuations and the understanding of their role are 
thus key points to develop a theory beyond MF, able to 



TABLE I: Critical points and critical exponents for the CP 
on UCM networks, compared with MF values. Numbers in 
parentheses indicate the uncertainty in the last digit. 
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0.63(4) 
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FIG. 4: Ratio of the standard deviation Apk of the distribu- 
tion of the density of occupied vertices of degree k over the 
average pk, versus the degree k, for 7 = 2.75, p = p c — 0.4215 
and two values of the system size N. 

fully capture the physics of collective phenomena on SF 
networks. 
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